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This  report  was  prepared  for  the  Harry  Diamond  Laboratories,  Adelphi , 
Maryland,  by  The  RDM  Corporation  under  Contract  DAAG39-75-C-0158.  This 
document  constitutes  the  final  teclinical  report  on  two  tasks  related  to 
improvements  in  the  NLT-2  Network  Analysis  Program. 


1. 


NONLINl'IAR  RESISTOR  MODEL 


A model  for  a nonlinear  resistor  has  been  developed  and  implemented  in 
NET-2.  The  nonlinear  resistor  has  actually  been  implemented  as  an 
algebraic  power  series  in  which  element  voltage  is  a function  of  element 
current.  The  power  series  function  is  defined  for  non-negative  values 
of  current  and  reflected  into  the  other  lialf  plane  (non-positive  current) 
as  a negative  function. 

C 

Details  of  the  power  series  element  are  given^on  the  following  page.  The 
FORTRAN  listing  of  subroutine  PSVI,  the  power  series  element  routine,  is 
given  in  Appendix  A. 


Power  Series  V-I  Curve 


Tlie  PSVI  element  is  a ircuit  element  which  may  be  used  to  represent  a 
voltage-current  characteristic  curve  in  which  the  voltage  across  the 
element  is  represented  by  a power  series  expansion  of  the  current  flowing 
tlirough  the  element.  The  voltage  across  the  element  is  related  to  the 
current  flowing  through  the  element  by  the  expression; 


V = 0 E K J 1 1 i 

i=l 

where:  V = element  voltage 

I = element  current 

^1’  ^2’  •••  *^n  ~ power  series  coefficients 
0 = sign  of  I 

From  the  above  equation  it  is  seen  that  the  element  voltage  for  negative 
values  of  current  is  the  negative  of  the  element  voltage  for  positive 
values  of  current,  i.e.,  if  V = F(l)  then  V = -F(-I)  for  I > 0.  Thi  s means 
that  the  PSVI  element  behaves  exactly  like  a nonlinear  resistor  whose  resis- 
tance is  given  by: 


r = ^=E 

i=l 


The  format  for  the  element  is: 

PSVIn  (p)  a b K,  , K , K 

12  n 

where:  PSVIn  = element  ID 

p = optional  parallel  parts  designation 
a and  b = node  pair  names 

Ki , , •••.  “ power  series  coefficients 

Ibe  power  series  coefficients  may  be  represented  by  numbers  or  mathematical 
expressions.  The  element  does  not  contain  any  intrinsic  computational  delay. 


The  voltage  across  the  element  is  given  by 

V = e - e, 
a b 

where  and  e-^j  are  the  voltages  at  nodes  a and  b,  respectively.  The  current 
flow  through  the  element  is  given  by 

I = V/R 

where  R is  the  equivalent  nonlinear  resistance  of  the  element,  element 
power  dissipation  P is  given  by 

P = IV 

The  user  may  I'eference  the  power  series  coefficients  , K2,  K^,  etc.,  by 
using  the  respective  symbolic  names  PSVIn,  PSVIn.l,  PSVIn.2,  etc.  The  element 
voltage,  current  and  power  dissipation  may  be  referenced  by  using  the  symbolic 
names  V(PSVIn),  l(PSVIn),  and  P (PSVIn),  respectively. 
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2.  OPTIMim'ION  IMPR0VI'MF:UTS 


Changes  have  been  made  to  the  optimization  routine  in  NET-2  to  permit  the 
user  to  obtain  output  information  from  a partially  completed  optimization 
solution  when  a computer  time  limit  has  been  reached.  The  results  printed 
are  the  best  optimized  values  which  have  been  found  up  to  the  moment  when 
the  time  limit  is  reached.  These  output  values  include  the  objective 
function  value,  the  values  of  the  free  par;ameters , and  the  values  of 
parameters  specified  by  constraint  relations. 

This  output  information  is  identical  to  that  which  is  provided  upon  normal 
completion  of  the  optimization  solution.  Ibe  optimization  results  are 
provided  whether  or  not  a checkpoint/restart  record  is  written  by  NET-2. 
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FORTRAN  Listing  of  Subroutine  PSVI 


SUBROUTINE  PSVI 
CUKMON/aU/Q(  1 ) 

CU^‘MDN/P/LnC(120),lXP(B97) 

CUf-‘lON/C277/L,LD,N,  I,  1X7700(3) 
CUhmDN/C277R/P»KJZ»V,RJ,D,S,RM,B,X7700(22) 
CUmmON/C278/1ALPHA,V1,IBETA 
Cn^M0N/C2bD/IRUN»X6D00(17) 

CU-M0N/C2  7 3/Y, y 1 , ICP, J, K,  1 X7301 
CU-moN/C200I/IX200,LAM,LAM1,LDYN,  1X201  (fa) 

C0MMDN/C2Lnc/LSAM,lXL0Cl(2),LDCN,lXL0C2(i(),L0CV,IXL0C3(2),L0Cb, 
1 lXLDt^(20) 

Cfil^MUN/C2faF/M,  IXfaFOO  ( 12) 

CL1»^MON/C2bDOl/IXfaDlO(fa),lCALC»IXbDll(3) 

COHMUN/C2b£/ICONV,IX6EOO(2) 

CU!'' MON/C  1 0800/XT  wn(  9 ),TwUlfa»XTvs02 
D I HE  MS  ION  ITAP 

data  ITAB/1 , ' 

EQUIVALENCE  LOC(llB)), (J0119,L0C(119)) 

CALL  ENTRYF  u 

*if(nut  res  _ variable  CALCULATIGN)THEN 
1F(IRUN,£Q,8)GU  TO  8000 

*DO(FDR  ALL  PSVI  ELEMENTS) 

L = Q(  J0U8  + LAM) 

LsL+JGl 19 

*INITIAL  SETUP 
LDsLDYN+TW01bA(0(L)-118,) 

JsQlL+l  ) 

J=0(LOCN+J) 

K=0(L+2) 

KrQ(LOCNf K) 

N=0(Lt3) 

aCASE(IRUN) 

GO  TU(1000» 9000, 3000»9000»5000ffa000»9000» 9000, b500)» IRON 
AlRUNsl,  STORE  B and  H TERMS 
P=l, 

IF(AND(0(L),10B).NE,0)P=a(L  + N + <4) 

RJZsARS(Q(LD) ) 

RM=N*0(L+N+3) 

IF  (N.f  Q,  1 )(;□  TO  laoo 
N = N-1 

RM  = RMARJZ  + ,M*G)(L  + Nt3) 

GO  TO  1300 
RM=ABS(RM)t2, 

YsP/KM 

•STORE  H TERMS 
CALL  STORY 
U(LDt2)=RM 


A-] 


c *1F(TIME  DOMAIN)ThEN 

1F(ICALC.NE,0)GU  TO  7000 
C ■‘STOKE  B TEKHS 

B=Y*Q(LDtl )-P*aCLD) 

IF(Q(LGAM  + J)  ,GT.O.  )(J(LOC3  + J)=0(LOCB  + J)tB 
IF(Q(LGAM  + K).GT.0,)O(LOCB-fK)  = Q(LDCB  + K)-B 
GO  TO  7000 

C *IHUN=i,  INITIALIZE  DC 

3000  G(LD)=0, 

Q(LD+1 )=0, 

GO  TO  7000 

C *lKUN=b,  CONVERGENCE  CHECK 

SOOO  Vl=faLUCVtJ)-0(LOCV  + K) 

KJ  = Q(l.D)  + (Vl-0(LD+l))/Q(LDt2I 
C *CALCULATE  V = F ( J) 

v = 0 

S=SI GN( 1 , , KJ) 

KjZrAHS(RJ) 

I = N 

SlOO  IF ( I .EG.O)GO  TO  5200 

v=  ( V>.J(L  + I + 3)  ) *KJZ 

1 = 1-1 

GO  TO  5100 

5200  v=S*V 

IF ( ABS( VI ) ,LT , lE-1 0, AND, A3S( V) ,LT , 1£-10)GO  TO  5«00 
D=ABS(V1+V) 

IF  CD.EQ, 0)GU  TO  5300 

IF(A0S(V1-V)/D.LE.  .OODGO  TO  5a00 
5300  1CONV=0 

5«00  aCLD)=KJ 

a(LDtn  = v 

GO  TO  7000 

C *1RUN=6,  sparse  matrix  setup 

6000  CALL  SKELl ( ITABC 1 ) ,2) 

GO  TO  7000 

C ■‘IRUNsR,  nodal  CONNECTIVITY 

6500  IF(M,NE.0)RETURN 

call  NODCDN(lTAB(n,2,0) 

GO  TO  7000 

7000  IF (LAM.EO.LAMl )GO  TO  9000 

cam=lam+i 
GO  TO  500 

C ^ELSE 

C ‘CALCULATE  RESPONSE  VARIABLES 

8000  V1=Q(LD+1) 

IF(IALPHA.EU.10)GO  TO  9000 

IF(IALPHA,EQ,12)V1=V1‘0(LD) 

IFdALPHA.EQ.  1 1)V1=Q(LD) 

9000  RETURN 
END 
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